axis.text.y = element_text(margin=margin(3,5,0,3,"pt")),
axis.ticks.length = unit(-0.7, "mm"),
text=element_text(size = text.size, family="Times"))
# export
####################################################
setwd(dir)
gA = ggplotGrob(beta_firm)
gB = ggplotGrob(beta_slave)
pdf("beta_merge.pdf", width = 7, height = 6)
grid.arrange( rbind(gA, gB, size = "last") )
dev.off()
beta_firm = ggplot() +
geom_vline(xintercept = conf_interval[1], linetype = 2, col = "grey50") +
geom_vline(xintercept = conf_interval[2], linetype = 2, col = "grey50") +
geom_density(data = beta_boot, aes(V1, fill = Source), alpha = 0.3, adjust = 2) +
scale_x_continuous(breaks = seq(0, 1.1, 0.1)) +
labs(x = expression(beta), y = "Density") +
coord_cartesian(xlim = c(0.15, 1.1), ylim = c(0, 60)) +
ggtitle("B.  United States Slave Estate") +
theme_bw() +
theme(panel.border = element_rect(color = "black"),
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
plot.title = element_text(face="bold", size = rel(1), hjust = 0.5),
legend.key.size = unit(0.4, "cm"),
axis.line = element_line(color = "black"),
axis.title.x = element_text(size = rel(0.9)),
axis.title.y = element_text(vjust= 1.1, size=rel(0.9)),
axis.text.x = element_text(margin=margin(5,5,0,0,"pt"), size = rel(1.2)),
axis.text.y = element_text(margin=margin(3,5,0,3,"pt")),
axis.ticks.length = unit(-0.7, "mm"),
text=element_text(size = text.size, family="Times"))
beta_slave = ggplot() +
geom_density(data = d, aes(beta, fill = Method), adjust = 1.5, alpha = 0.4) +
scale_x_continuous(breaks = seq(0, 1.1, 0.1)) +
labs(x = expression(beta), y = "Density") +
coord_cartesian(xlim = c(0.15, 1.1), ylim = c(0, 10)) +
ggtitle("A.  Case Study Firms") +
theme_bw() +
theme(panel.border = element_rect(color = "black"),
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
plot.title = element_text(face="bold", size = rel(1), hjust = 0.5),
legend.key.size = unit(0.4, "cm"),
axis.line = element_line(color = "black"),
axis.title.x = element_text(size = rel(0.9)),
axis.title.y = element_text(vjust= 1.1, size=rel(0.9)),
axis.text.x = element_text(margin=margin(5,5,0,0,"pt"), size = rel(1.2)),
axis.text.y = element_text(margin=margin(3,5,0,3,"pt")),
axis.ticks.length = unit(-0.7, "mm"),
text=element_text(size = text.size, family="Times"))
beta_firm = ggplot() +
geom_vline(xintercept = conf_interval[1], linetype = 2, col = "grey50") +
geom_vline(xintercept = conf_interval[2], linetype = 2, col = "grey50") +
geom_density(data = beta_boot, aes(V1, fill = Source), alpha = 0.3, adjust = 2) +
scale_x_continuous(breaks = seq(0, 1.1, 0.1)) +
labs(x = expression(beta), y = "Density") +
coord_cartesian(xlim = c(0.15, 1.1), ylim = c(0, 60)) +
ggtitle("B.  United States Slave Estate") +
theme_bw() +
theme(panel.border = element_rect(color = "black"),
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
plot.title = element_text(face="bold", size = rel(1), hjust = 0.5),
legend.key.size = unit(0.4, "cm"),
axis.line = element_line(color = "black"),
axis.title.x = element_text(size = rel(0.9)),
axis.title.y = element_text(vjust= 1.1, size=rel(0.9)),
axis.text.x = element_text(margin=margin(5,5,0,0,"pt"), size = rel(1.2)),
axis.text.y = element_text(margin=margin(3,5,0,3,"pt")),
axis.ticks.length = unit(-0.7, "mm"),
text=element_text(size = text.size, family="Times"))
# export
####################################################
setwd(dir)
gA = ggplotGrob(beta_firm)
gB = ggplotGrob(beta_slave)
pdf("beta_merge.pdf", width = 7, height = 6)
grid.arrange( rbind(gA, gB, size = "last") )
dev.off()
beta_firm = ggplot() +
geom_vline(xintercept = conf_interval[1], linetype = 2, col = "grey50") +
geom_vline(xintercept = conf_interval[2], linetype = 2, col = "grey50") +
geom_density(data = beta_boot, aes(V1, fill = Source), alpha = 0.3, adjust = 2) +
scale_x_continuous(breaks = seq(0, 1.1, 0.1)) +
labs(x = expression(beta), y = "Density") +
coord_cartesian(xlim = c(0.15, 1.1), ylim = c(0, 60)) +
ggtitle("A.  United States Slave Estate") +
theme_bw() +
theme(panel.border = element_rect(color = "black"),
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
plot.title = element_text(face="bold", size = rel(1), hjust = 0.5),
legend.key.size = unit(0.4, "cm"),
axis.line = element_line(color = "black"),
axis.title.x = element_text(size = rel(0.9)),
axis.title.y = element_text(vjust= 1.1, size=rel(0.9)),
axis.text.x = element_text(margin=margin(5,5,0,0,"pt"), size = rel(1.2)),
axis.text.y = element_text(margin=margin(3,5,0,3,"pt")),
axis.ticks.length = unit(-0.7, "mm"),
text=element_text(size = text.size, family="Times"))
beta_slave = ggplot() +
geom_density(data = d, aes(beta, fill = Method), adjust = 1.5, alpha = 0.4) +
scale_x_continuous(breaks = seq(0, 1.1, 0.1)) +
labs(x = expression(beta), y = "Density") +
coord_cartesian(xlim = c(0.15, 1.1), ylim = c(0, 10)) +
ggtitle("A.  Antebellum US Slave Estate") +
theme_bw() +
theme(panel.border = element_rect(color = "black"),
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
plot.title = element_text(face="bold", size = rel(1), hjust = 0.5),
legend.key.size = unit(0.4, "cm"),
axis.line = element_line(color = "black"),
axis.title.x = element_text(size = rel(0.9)),
axis.title.y = element_text(vjust= 1.1, size=rel(0.9)),
axis.text.x = element_text(margin=margin(5,5,0,0,"pt"), size = rel(1.2)),
axis.text.y = element_text(margin=margin(3,5,0,3,"pt")),
axis.ticks.length = unit(-0.7, "mm"),
text=element_text(size = text.size, family="Times"))
beta_firm = ggplot() +
geom_vline(xintercept = conf_interval[1], linetype = 2, col = "grey50") +
geom_vline(xintercept = conf_interval[2], linetype = 2, col = "grey50") +
geom_density(data = beta_boot, aes(V1, fill = Source), alpha = 0.3, adjust = 2) +
scale_x_continuous(breaks = seq(0, 1.1, 0.1)) +
labs(x = expression(beta), y = "Density") +
coord_cartesian(xlim = c(0.15, 1.1), ylim = c(0, 60)) +
ggtitle("A.  Case-Study Firm") +
theme_bw() +
theme(panel.border = element_rect(color = "black"),
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
plot.title = element_text(face="bold", size = rel(1), hjust = 0.5),
legend.key.size = unit(0.4, "cm"),
axis.line = element_line(color = "black"),
axis.title.x = element_text(size = rel(0.9)),
axis.title.y = element_text(vjust= 1.1, size=rel(0.9)),
axis.text.x = element_text(margin=margin(5,5,0,0,"pt"), size = rel(1.2)),
axis.text.y = element_text(margin=margin(3,5,0,3,"pt")),
axis.ticks.length = unit(-0.7, "mm"),
text=element_text(size = text.size, family="Times"))
# export
####################################################
setwd(dir)
gA = ggplotGrob(beta_firm)
gB = ggplotGrob(beta_slave)
pdf("beta_merge.pdf", width = 7, height = 6)
grid.arrange( rbind(gA, gB, size = "last") )
dev.off()
beta_slave = ggplot() +
geom_density(data = d, aes(beta, fill = Method), adjust = 1.5, alpha = 0.4) +
scale_x_continuous(breaks = seq(0, 1.1, 0.1)) +
labs(x = expression(beta), y = "Density") +
coord_cartesian(xlim = c(0.15, 1.1), ylim = c(0, 10)) +
ggtitle("B.  Antebellum US Slave Estate") +
theme_bw() +
theme(panel.border = element_rect(color = "black"),
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
plot.title = element_text(face="bold", size = rel(1), hjust = 0.5),
legend.key.size = unit(0.4, "cm"),
axis.line = element_line(color = "black"),
axis.title.x = element_text(size = rel(0.9)),
axis.title.y = element_text(vjust= 1.1, size=rel(0.9)),
axis.text.x = element_text(margin=margin(5,5,0,0,"pt"), size = rel(1.2)),
axis.text.y = element_text(margin=margin(3,5,0,3,"pt")),
axis.ticks.length = unit(-0.7, "mm"),
text=element_text(size = text.size, family="Times"))
# export
####################################################
setwd(dir)
gA = ggplotGrob(beta_firm)
gB = ggplotGrob(beta_slave)
pdf("beta_merge.pdf", width = 7, height = 6)
grid.arrange( rbind(gA, gB, size = "last") )
dev.off()
# export
####################################################
setwd(dir)
gA = ggplotGrob(beta_firm)
gB = ggplotGrob(beta_slave)
pdf("beta_merge.pdf", width = 6, height = 6)
grid.arrange( rbind(gA, gB, size = "last") )
dev.off()
# export
####################################################
setwd(dir)
gA = ggplotGrob(beta_firm)
gB = ggplotGrob(beta_slave)
pdf("beta_merge.pdf", width = 7, height = 6)
grid.arrange( rbind(gA, gB, size = "last") )
dev.off()
beta_slave = ggplot() +
geom_density(data = d, aes(beta, fill = Method), adjust = 1.5, alpha = 0.4) +
scale_x_continuous(breaks = seq(0, 1.1, 0.1)) +
labs(x = expression(beta), y = "Density") +
coord_cartesian(xlim = c(0.15, 1.1), ylim = c(0, 10)) +
ggtitle("B.  Antebellum US Slave Estate") +
theme_bw() +
theme(panel.border = element_rect(color = "black"),
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
plot.title = element_text(face="bold", size = rel(1), hjust = 0.5),
legend.key.size = unit(0.4, "cm"),
legend.position = c(0.8, 0.7),
axis.line = element_line(color = "black"),
axis.title.x = element_text(size = rel(0.9)),
axis.title.y = element_text(vjust= 1.1, size=rel(0.9)),
axis.text.x = element_text(margin=margin(5,5,0,0,"pt"), size = rel(1.2)),
axis.text.y = element_text(margin=margin(3,5,0,3,"pt")),
axis.ticks.length = unit(-0.7, "mm"),
text=element_text(size = text.size, family="Times"))
beta_firm = ggplot() +
#geom_vline(xintercept = conf_interval[1], linetype = 2, col = "grey50") +
#geom_vline(xintercept = conf_interval[2], linetype = 2, col = "grey50") +
geom_density(data = beta_boot, aes(V1, fill = Source), alpha = 0.3, adjust = 2) +
scale_x_continuous(breaks = seq(0, 1.1, 0.1)) +
labs(x = expression(beta), y = "Density") +
coord_cartesian(xlim = c(0.15, 1.1), ylim = c(0, 60)) +
ggtitle("A.  Case-Study Firm") +
theme_bw() +
theme(panel.border = element_rect(color = "black"),
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
plot.title = element_text(face="bold", size = rel(1), hjust = 0.5),
legend.key.size = unit(0.4, "cm"),
legend.position = c(0.8, 0.5),
axis.line = element_line(color = "black"),
axis.title.x = element_text(size = rel(0.9)),
axis.title.y = element_text(vjust= 1.1, size=rel(0.9)),
axis.text.x = element_text(margin=margin(5,5,0,0,"pt"), size = rel(1.2)),
axis.text.y = element_text(margin=margin(3,5,0,3,"pt")),
axis.ticks.length = unit(-0.7, "mm"),
text=element_text(size = text.size, family="Times"))
# export
####################################################
setwd(dir)
gA = ggplotGrob(beta_firm)
gB = ggplotGrob(beta_slave)
pdf("beta_merge.pdf", width = 7, height = 6)
grid.arrange( rbind(gA, gB, size = "last") )
dev.off()
gA = ggplotGrob(beta_firm)
gB = ggplotGrob(beta_slave)
pdf("beta_merge.pdf", width = 6, height = 6)
grid.arrange( rbind(gA, gB, size = "last") )
dev.off()
beta_firm = ggplot() +
#geom_vline(xintercept = conf_interval[1], linetype = 2, col = "grey50") +
#geom_vline(xintercept = conf_interval[2], linetype = 2, col = "grey50") +
geom_density(data = beta_boot, aes(V1, fill = Source), alpha = 0.3, adjust = 2) +
scale_x_continuous(breaks = seq(0, 1.1, 0.1)) +
labs(x = expression(beta), y = "Density") +
coord_cartesian(xlim = c(0.15, 1.1), ylim = c(0, 60)) +
ggtitle("A.  Case-Study Firm") +
theme_bw() +
theme(panel.border = element_rect(color = "black"),
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
plot.title = element_text(face="bold", size = rel(1), hjust = 0.5),
legend.key.size = unit(0.4, "cm"),
legend.position = c(0.8, 0.5),
axis.line = element_line(color = "black"),
axis.title.x = element_text(size = rel(0.9)),
axis.title.y = element_text(vjust= 1.1, size=rel(0.9)),
axis.text.x = element_text(margin=margin(5,5,0,0,"pt"), size = rel(1.2)),
axis.text.y = element_text(margin=margin(3,5,0,3,"pt")),
axis.ticks.length = unit(-0.7, "mm"),
text=element_text(size = text.size, family="Times"))
# export
####################################################
setwd(dir)
gA = ggplotGrob(beta_firm)
gB = ggplotGrob(beta_slave)
pdf("beta_merge.pdf", width = 6, height = 6)
grid.arrange( rbind(gA, gB, size = "last") )
dev.off()
source('~/Desktop/origin_inequality/Supplementary Material/Figures/beta_slave.R')
# export
####################################################
setwd(dir)
gA = ggplotGrob(beta_firm)
gB = ggplotGrob(beta_slave)
pdf("beta_merge.pdf", width = 6, height = 6)
grid.arrange( rbind(gA, gB, size = "last") )
dev.off()
source('~/Desktop/origin_inequality/Supplementary Material/Figures/beta_slave.R')
# export
####################################################
setwd(dir)
gA = ggplotGrob(beta_firm)
gB = ggplotGrob(beta_slave)
pdf("beta_merge.pdf", width = 6, height = 6)
grid.arrange( rbind(gA, gB, size = "last") )
dev.off()
source('~/Desktop/origin_inequality/Supplementary Material/Figures/pay_concentration_gini_empirical.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/pay_concentration_gini_empirical.R')
gini = ggplot() +
#geom_point(data = mod, aes(x = energy, y = power_gini), size = 0.2, col = "grey30") +
geom_point(data = mod, aes(x = energy, y = pay_gini, color = pay_exponent), size = 0.2) +
geom_point(data = d, aes(x = EG.USE.PCAP.KG.OE, y = SI.POV.GINI)) +
#geom_polygon(data = frontier, aes(x = energy, y = top), fill = "black", alpha = 0.2, col = "black") +
geom_point(data = milanovic, aes(x = fit, y = Gini1/100)) +
geom_point(data = milanovic, aes(x = fit, y = Gini2/100)) +
scale_x_log10("Energy Use per Capita (GJ/year)",  breaks = c(2, 5, 10, 20, 50, 100, 200, 500, 1000, 2000)) +
scale_y_continuous("Gini Index", breaks = seq(0,1, 0.1)) +
scale_color_gradientn(expression(beta), colours = rainbow(8), breaks = seq(0, 1, 0.1), labels = seq(0, 1, 0.1)) +
coord_cartesian(xlim = seq(2.5, 2000), ylim = c(0, 1)) +
theme_bw() +
theme(panel.border = element_rect(color = "black"),
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
plot.title = element_text(face="bold", size = rel(1), hjust = 0.5),
legend.title = element_text(hjust = 0.2),
legend.key.height = unit(1.2, "cm"),
legend.justification = c(0, 0.3),
axis.line = element_line(color = "black"),
axis.title.x=element_text(vjust=-0.3, size=rel(0.9)),
axis.title.y=element_text(vjust= 1.1, size=rel(0.9)),
axis.text.x = element_text(margin=margin(5,5,0,0,"pt")),
axis.text.y = element_text(margin=margin(3,5,0,3,"pt")),
axis.ticks.length = unit(-0.7, "mm"),
text=element_text(size = text.size, family="Times"))
# export
####################################################
setwd(dir)
gA = ggplotGrob(gini)
png("pay_concentration_gini_empirical.png", width = 7, height = 4,  units = 'in', res = 600)
grid.arrange(gA)
dev.off()
source('~/Desktop/origin_inequality/Supplementary Material/Figures/pay_concentration_gini_empirical.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/pay_concentration_gini_empirical.R')
library(ggplot2)
library(gridExtra)
library(data.table)
library(hmod)
library(magrittr)
library(scales)
library(here)
library(ineq)
text.size = 10
# empirical
dir = here()
wd = gsub("Figures", "Empirical Data/World Bank", dir)
setwd(wd)
d = fread("data_format.csv") %>% na.omit()
d$SI.POV.GINI = d$SI.POV.GINI/100
d$EG.USE.PCAP.KG.OE = d$EG.USE.PCAP.KG.OE * 0.041868 + 3 # add food energy
# remove island states
remove = c("Comoros", "Kiribati", "Solomon Islands", "Cabo Verde")
d = d[! country %in% remove]
# inequality frontier
####################################################################################
top_share = function(energy_subsistence, energy_pc, top_n){
n_people = 100/top_n
energy_bottom = energy_subsistence*(n_people - 1)
energy_total = energy_pc*n_people
energy_top = energy_total - energy_bottom
population = c(rep(energy_subsistence, n_people-1), energy_top)
g = ineq(population)
return(g)
}
kcal_to_gj = 0.000004186798
top_n = 1
energy_subsistence_low = kcal_to_gj*2000*365
energy_low = exp(seq(log(energy_subsistence_low), log(3000), length.out = 1000))
top_low = mapply(top_share, energy_subsistence_low, energy_low, top_n)
energy_subsistence_high = kcal_to_gj*3000*365
energy_high = exp(seq(log(energy_subsistence_high), log(3000), length.out = 1000))
top_high = mapply(top_share, energy_subsistence_high, energy_high, top_n)
frontier = rbind(
data.table(energy = energy_low, top = top_low),
data.table(energy = energy_high, top = top_high)[order(-top)],
data.table(energy = energy_low, top = top_low)[1,]
)
# milanovic
#####################################
wd = gsub("Figures", "Empirical Data/Milanovic", dir)
setwd(wd)
milanovic = fread("milanovic_energy.csv")
# model
##################################################
dir = here()
wd = gsub("Figures", "Hierarchy Model/results", dir)
setwd(wd)
mod = fread("inequality.csv")
names(mod) = c("firm_size", "power_gini", "pay_gini", "pay_exponent", "power_top_1", "pay_top_1")
wd = gsub("Figures", "Empirical Data/GEM", dir)
setwd(wd)
energy = fread("energy_v_firm.csv")
x = log(energy$energy)
y = log(energy$firm_mean)
r = lm(y ~ x)
a = exp(coef(r)[1])
b = coef(r)[2]
e_predict = function(x){
y = (x/a)^(1/b)
return(y)
}
energy_predict = e_predict(mod$firm_size)
mod$energy = energy_predict
# historical energy data
kcal_to_gj = 0.000004186798
metabolic_low = kcal_to_gj*2000*365
metabolic_high = kcal_to_gj*2500*365
western_BCE_14000 = kcal_to_gj*365*4000 # Morris p. 61
western_BCE_10000 = kcal_to_gj*365*5000 # Morris p. 61
eastern_3000_bce = kcal_to_gj*365*8000  # Morris p. 111
eastern_1500_ce = kcal_to_gj*365*30000  # Morris p. 111
US_min = 133
US_max = 379
Qatar_min = 770
Qatar_max = 953
geo_mean = function(x,y){
exp(mean(log(c(x,y))))
}
View(milanovic)
source('~/Desktop/origin_inequality/Supplementary Material/Figures/pay_concentration_gini_empirical.R')
source('~/Desktop/origin_inequality/Supplementary Material/Figures/pay_concentration_gini_empirical.R')
library(ggplot2)
library(gridExtra)
library(data.table)
library(hmod)
library(magrittr)
library(scales)
library(here)
library(Hmisc)
text.size = 10
dir = here()
wd = gsub("Figures", "Hierarchy Model/results", dir)
setwd(wd)
mod = fread("inequality.csv")
names(mod) = c("firm_size", "power_gini", "pay_gini", "pay_exponent")
wd = gsub("Figures", "Empirical Data/GEM", dir)
setwd(wd)
energy = fread("energy_v_firm.csv")
x = log(energy$energy)
y = log(energy$firm_mean)
r = lm(y ~ x)
a = exp(coef(r)[1])
b = coef(r)[2]
e_predict = function(x){
y = (x/a)^(1/b)
return(y)
}
energy_predict = e_predict(mod$firm_size)
mod$energy = energy_predict
# historical energy data
kcal_to_gj = 0.000004186798
metabolic_low = kcal_to_gj*2000*365
metabolic_high = kcal_to_gj*2500*365
western_BCE_14000 = kcal_to_gj*365*4000 # Morris p. 61
western_BCE_10000 = kcal_to_gj*365*5000 # Morris p. 61
eastern_3000_bce = kcal_to_gj*365*8000  # Morris p. 111
eastern_1500_ce = kcal_to_gj*365*30000  # Morris p. 111
kcal_to_gj*365*25000
western_BCE_14000 = kcal_to_gj*365*4000
western_BCE_10000 = kcal_to_gj*365*5000
kcal_to_gj*365*4,438
source('~/.active-rstudio-document')
kcal_to_gj*365*4438
kcal_to_gj*365*15500
kcal_to_gj*365*43000
kcal_to_gj*365*5864
kcal_to_gj*365*10710
kcal_to_gj*365*10710
kcal_to_gj*365*20000
kcal_to_gj*365*25000
kcal_to_gj*365*25000
kcal_to_gj*365*4000
kcal_to_gj*365*30000
kcal_to_gj*365*27000
kcal_to_gj*365*2500
kcal_to_gj*365*24000
energy = fread("wdi.csv")
source('~/Desktop/origin_inequality/Supplementary Material/Empirical Data/GEM/GEM_energy_match.R')
energy = fread("wdi.csv")
wd = gsub("GEM", "World Bank", dir)
setwd(wd)
energy = fread("wdi.csv")
source('~/Desktop/origin_inequality/Supplementary Material/Empirical Data/GEM/GEM_energy_match.R')
energy = energy[is.na(energy_pc_gj) == F]
energy$year_country = paste(energy$country.code, energy$year)
energy = energy[is.na(energy_pc_gj) == F]
setwd(dir)
GEM = fread("GEM_mean_firm_size.csv")
GEM$year_country = paste(GEM$code, GEM$year)
GEM_energy_quantiles = fread("GEM_energy_quantiles.csv")
merge = merge(GEM, energy, by = "year_country", no.dups = T)
merge = merge[(order(country.x))]
GEM_energy_quantiles = GEM_energy_quantiles[country %in% merge$country.x]
GEM_energy_quantiles = GEM_energy_quantiles[order(country)]
merge$energy_quantile = GEM_energy_quantiles$energy_quantile
source('~/Desktop/origin_inequality/Supplementary Material/Empirical Data/GEM/GEM_energy_match.R')
merge = merge(GEM, energy, by = "year_country", no.dups = T)
View(merge)
View(energy)
energy$year_country = paste(code3, energy$year)
energy$year_country = paste(energy$code3, energy$year)
energy = energy[is.na(energy_pc_gj) == F]
source('~/Desktop/origin_inequality/Supplementary Material/Empirical Data/GEM/GEM_energy_match.R')
result = data.table(country = merge$country.x,
code = merge$code,
year = merge$year.x,
firm_lower = merge$lower,
firm_mean = merge$firm_mean,
firm_upper = merge$upper,
energy = merge$energy_pc_gj,
energy_quantile = merge$energy_quantile)
source('~/Desktop/origin_inequality/Supplementary Material/Empirical Data/GEM/GEM_energy_match.R')
View(merge)
View(result)
